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Abstract 

The probability to detect a neutrino produced in pion decay at a 
finite distance exhibits unique interference properties that depend on 
the absolute mass of the neutrino. We describe the neutrino, lepton, 
and pion by a many-body wave function and find the probability is 
subject to a large finite-size correction. Its rate at a distance L is 
expressed as Tq + g(u} u L/c)Ti, where g(u u L/c) is the universal func- 
tion, lo v = m1c A / '(2E u h), c is the speed of light, and Fq is a constant 
computed with the standard plane-wave S-matrix. The finite-size cor- 
rection is rigorously computed using wave packets via the light-cone 
singularity of a system composed of the pion and charged lepton and 
reveals the diffraction pattern of a single quantum. We discuss the 
implications of this correction for the muon-neutrino and electron- 
neutrino reactions. With sufficient statistics, the neutrino diffraction 
would supply the absolute mass of the neutrino. 
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1 Neutrino interference 



Interference phenomena involving photons, electrons, neutrons, and other 
heavy elements are useful for testing quantum mechanics and other basic 
principles. Neutrinos are elementary particles that are produced in weak de- 
cays of pions and other weak processes and are described by quantum waves 
that follow the superposition principle. Since neutrinos interact extremely 
weakly with matter, the probability to detect neutrinos at a finite distance 
exhibits a single quantum interference that makes the probability receive 
an unusual finite-size correction. This pattern could, via its unique behav- 
ior, provide an absolute value of the neutrino mass, which is at present an 
unknown physical quantity. 

A system composed of a pion, charged lepton, and neutrino are described 
by the following Lagrangian of a field of pion, <p(x), of charged lepton, l(x), 
and of neutrino, u(x), 

L = Lq + Li n t, 

Lq = d^d^ip - mlip*ip + l(j-p- mi)l + u(j ■ p - m v )v, 

Lint = gJL'aL x JIZL 9 = G F /V2, (1) 



where Gf is the Fermi coupling constant and JY~ A is V — A current. The 
system is invariant under a Lorentz transformation and momenta can be any 
value from — oo to +oo. A pion and decay products are expressed by a many 
body wave function |\&) that follows a Schrodinger equation 

ih^\*) = (H + H int )\V), (2) 

where the free part, H , and the interaction part, H int , are derived from the 
previous Lagrangian. A time-dependent solution in the first order of H int is 

l*> = e"*** (j^ (0) > + J d(5D(u,t)\(5)((5\H mt \^) 
u = Ep- E , H \f3) = Ep\P), H \^) = E \i/,W), (3) 
where d/3 is a measure for a complete set of \(3), and 

e^ir* - 1 

D(u,t) = . (4) 

oo 
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In Eq. 02]), all the possible states of \f3) that couple with H int must be in- 
cluded. 

D(oj,t) at t = oo and an average over a finite interval, At of uAt ^> 1, 
behave as 

D(u,t) = -2txi8(uj)\ t ->■ oo, (5) 
average D(u,t) = ; finite £. (6) 

Thus the wave function at t = oo and the average over a finite interval at a 
finite t are 

I**) = e-**« (V (0) > - 2m J d(36(E p - E )\P)(p\H int \^^J , (7) 

^average) = ^ (|^ (0) ) - | \0) (/3|^ nt |^ (0) )) , (8) 

and satisfy in the first order of i? ini , 

tfltfoo) = £ |^oo>, ifol^oc) = ^Ol^oo), (9) 
H I ^average) Eq | ^average) , -^0 1 ^ average) 7^ -^0 | ^a-uerage) • (l 0) 

In both cases, the frequency is E /h and the total energy defined by H is E . 
The total energy is conserved. The state of a momentum eigenstate at t = 

l^ (0) ) = (11) 

evolves according to Eq. fl3]) and is a coherent superposition of and 
plane waves of a lepton and neutrino with the particular weight. The state 
at a finite t is different from that at t — oo, which is used for an ordinary 
calculation of the decay rate. 

At the infinite time t, the state becomes a superposition of \ip^) and 
| /So) °f the same kinetic energy Ep Q = E . Since the state |/3q) has the same 
kinetic energy, this is regarded as a state composed of free particles. The 
decay products follow dynamics of classical particles, and their norm shows 
the probability, which increases linearly with t, hence a pion decays with a 
constant rate like a classical particle. 

At a finite time t, on the other hand, the state is a superposition of 
and those of an arbitrary value of Ep > 0. Since the kinetic energy Ep is 
different from E and not a constant, the states in this region are different 
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from the state composed of free particles. The particle picture does not hold, 
and the wave nature remains here. 

The evolution of the state Eq. fl3]) in time is similar to an evolution of a 
stationary state of a system with potential or obstacle in x, where the kinetic 
energy varies due to a potential energy. A superposition of plane waves is 
formed and reveals interference of waves of many components, which we call 
diffraction in this paper. A superposition of waves Eq. at a finite t could 
show a similar phenomenon. In fact, the kinetic energy of waves composed of 
a charged lepton and neutrino in Eq. (J3J) at a finite t is not constant. Hence 
this could show a non-uniform space-time dependent behavior, which is a 
feature of diffraction. Now the weight of states is determined by H int and 
this diffraction phenomenon should have a universal property. After passing 
this period, the state has the energy Eq and becomes to be particle-like again 
at t — > oo, in the same way as the stationary wave at |af| — >• oo. 

Physical quantities are measured by observing the lepton or neutrino. 
Probability to detect the particle is finite experimentally and a theoretical 
value must be finite and unique with an amplitude that satisfies boundary 
conditions of experiments. The neutrino is identified from its reactions with 
nucleus that are caused by the weak interaction and its properties are deter- 
mined by the probabilities that these processes occur. Hence the final state 
of the transition processes is expressed by a wave packet of the nuclear size. 
Accordingly the probability to detect this neutrino is expressed by the wave 
packets. Wave packets that satisfy free wave equations and are localized in 
space are important for rigorously defining scattering amplitude [II [2] • 

A standard method to calculate the decay rate that uses the transition 
amplitudes of plane waves with asymptotic boundary conditions at time 
t = ±00 gives the asymptotic quantities at T = 00, where T is a time 
interval between in and out states. But that gives neither quantities at a 
finite T nor a position-dependent probability. This is because the boundary 
conditions are different and the plane waves are translationally invariant. 
To compute the position-dependent probability, the wave function express- 
ing entire pion-decay process, which begins as a superposition of a one-pion 
state and the state of the decay products, and satisfying the boundary con- 
ditions are necessary. The amplitude of this situation is expressed by S[T] 
satisfying the boundary conditions at a finite T given below. 

S[T] is defined by M0ller operators at a finite T, fi±(T), as S[T] = 
fiL(T)fi + (T). f2±(T) are expressed by a free Hamiltonian H and a total 
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Hamiltonian H by fi±(T) = lim e lHt e °*. From this expression, S[T] 
satisfies 

[^[T],iyo]=z{^l(T)}o + (T)-zQl(T)^0 + (T). (12) 

Thus a kinetic energy defined by H is not conserved except at T = oo 
and a matrix element of S[T] between two eigenstates of H Q , \a) and \(3) of 
eigenvalues E a and Ep, is decomposed into two components 

(P\S[r\\a) = (P\S^[T)\a) + (^\S^[T]\a), (13) 

where satisfies Ep = E a and satisfies Ep ^ E a . At T — >■ oo, the 
right-hand side of Eq. ffl2|) vanishes and so does the second term of Eq. ffl3|) . 
Hence the energy defined by Hq is conserved. At a finite T, on the other 
hand, the right-hand side of Eq. ffl2|) and second term of Eq. fflBl do not 
vanish. The states of Ep ^ E a couple and give the finite-size correction. The 
states of the different eigenvalues of Hq are orthogonal each other and the 
cross term of the first and second terms of Eq. (1131) in a square of the modulus 
vanishes. Consequently the finite-size correction to the probability becomes 
positive semi-definite. Unitarity S'fTjS'^fT] = 1 is satisfied and ensures the 
conservation of probability. 

Innumerable states in the ultra-violet energy region in Eq. ([3]) propagate 
with the speed of light and couple to S[T], hence they necessarily give a 
finite-size correction to a particle with small mass. Although the finite-time- 
interval effect has been considered insignificant and ignored heretofore, the 
correction is actually non-negligible and important especially to recent and 
future neutrino experiments. We show herein that the probability to detect a 
neutrino is studied with S[T] and is affected by the large finite-size correction 
in a macroscopic area. The diffraction of a neutrino is a main part of the 
finite-size correction for a system of a pion and charged lepton, and neutrino 
wave, and becomes visible in macroscopic scales. 

Because an asymptotic boundary condition in scattering or decay pro- 
cesses [H [2] is satisfied manifestly with wave packets, textbooks on scattering 
theory [SI HI [5] emphasize the importance of wave packets. They give further- 
more proofs that the asymptotic values are computed with plane waves. For 
computing non-asymptotic value that depends on a time interval, the wave 
packets are necessary. In the process of detecting neutrino, the neutrino wave 
packet [HI QUI EEE] is well localized and expresses a nucleon wave function in 
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v(p u , X u , TV) 
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Figure 1: A pion of a momentum p n decays to a muon of a momentum p^ and 
a neutrino of a wave packets (p v , X„). Since the neutrino wave packet becomes 
finite along a narrow region of a velocity v u , a coordinate x is integrated over 
this region. This quantum mechanical amplitude resembles diffraction of a 
classical wave that passes through a hole. 

the nucleus with which the neutrino interacts [T^[T^[15lfT6lfT7tfT8|[T9tl20]. 
From the probability obtained with this amplitude, the physical quantities 
measured in experiments are calculated, such as the number of events, the 
flux, and others. The initial pion is expressed by a large wave packet, so the 
plane-wave approximation is very good and is therefore used in this paper. 
The unobserved muon is also expressed by a plane wave. 

Mass-squared differences 8ml are negligibly small [6], and their central 
value is currently unknown, so we study the situation in which the average- 
mass-squared ml satisfies, ml ^> bm 2 v . To begin, we present the one-flavor 
case. Extension to the general case is straightforward. 



We now compute the probabilities to detect neutrinos and charged leptons, 
I, at a finite distance. An amplitude to detect a neutrino of an average 
momentum p v at X„ in the decay of a pion prepared at t = T^ of a momentum 
p,,-, Fig. (JT]), is expressed as T = f d A x (I, u\H w (x)\tt) , where a lepton / has a 
momentum pi and a neutrino is expressed by a wave packet. These states are 
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Position-dependent probability 
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expressed as \n) = Ip^T^), |/, v) = z/, p v , X„, T u ). T is written with 

the hadronic matrix element of V — A current and Dirac spinors 



T = j (PxdkvNiiOlJ^MWuipfail-jsMk 
x exp 



-i{PTT ~Pi) ■ x/h + ik v ■ (x - X u )/h- -Vvf , (11) 



where Aq = igioyjix^ {mim u / EiE v )^ . The time t is integrated over the 
region TV < t < T u . a v is the size of the neutrino wave packet and is 
estimated from the size of a nucleus. For the sake of simplicity, we use the 
Gaussian form of the wave packet in this paper. The result for the finite-size 
correction is the same in general wave packets, as will be verified in Appendix 
C. The amplitude T depends on the momenta of the final state and also on 
coordinates (T^Xj,) and T^. 

Integrating over k u in Eq. (THj) . we obtain a Gaussian function of x about 
the moving center with the velocity v v = p u c 2 /E u , x = X u +v u (t—T u ), which 
vanishes at large \x — x \ and satisfies the asymptotic boundary condition. 
The integrand becomes finite along a narrow space-time region of the velocity 
v u , and is similar to classical waves that pass through a small hole, which 
show diffraction due to boundary conditions. If we integrate over the space 
coordinate x and the time t in the finite-time interval T = T u — T^, we find 
the amplitude proportional to 

s ^Ml exp f-^6p 2 ) , u = 6E-v„-6p, 
u V 2 / 

5E = Eiip^) + E v (p u ) - E 7T (p 7T ), Sp = p u +pi-p n . (15) 

In the limit a v and T — > oo, the amplitude is proportional to the four di- 
mensional delta function (2ir) 4 '5^(5p), and becomes non-zero in the finite 
kinematical region of the final states. If a v is finite, bp is not strictly zero 
and u deviates from 6E. At a large T, the amplitude gets dominant contri- 
butions from to w 0. Since the frequency u is 5E — v ■ 5p, which is equivalent 
to that of a system moving with velocity v u , u — has a solution of 5p ^ 
and 5E ^ (see Appendix A) in addition to the normal solution, 5p = and 
SE = 0. Thus the number of states at u w becomes much larger than that 
of SE w and the probability integrated over the final states approaches 
the asymptotic value very slowly with T. Hence the quantum mechanical 
amplitude to detect a neutrino at a finite T has a large finite-size correction. 
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The amplitude in quantum mechanics that includes the detection process, 
Eq. (j!4j) . is in fact almost equivalent to the classical waves of the diffraction 
caused by the boundary conditions. Accordingly it would be reasonable to 
call the present phenomenon "neutrino diffraction" . Being a purely quantum 
mechanical phenomenon, neutrino diffraction has various unique properties 
that are different from the diffraction of the classical waves. 

The total probability is an integral of a square of the modulus of the 
amplitude over the complete set of final states [9] 



S1,S2 



where the momenta of the neutrino and muon are integrated over whole 
positive energy region and the position of the wave packet is integrated over 
the region of the detector. This depends on a time interval T Hereafter the 
natural unit, c = h = 1, is taken in majority of places, but c and H are 
written explicitly when it is necessary. 

To rigorously compute the finite-size correction in a consistent manner 
with the Lorentz invariance, we express the probability Eq. (fT6|) with a cor- 
relation function A 7r ;(xi,x 2 ) of the coordinates. The total probability is 
finite and the order of integrations of the coordinates and momentum is in- 
terchangeable. Here the muon momentum is integrated first for fixed Xi, after 
summing over the spin, the probability is written in the form 



P= I d ^^y,~^ I d A x 1 d 4 x 2 exp 



— - x f) 2 + i<j>{8x) 



2a v 

x f = X u + v u (U - T u ), 5x = xi - x 2 , <j)(8x) = p u -5x, (17) 

3 

where N 2 = g 2 (Air / a u )^ V' 1 , V is a normalization volume for the initial 
pion, and 

A^Ofc) = JL | J|L { 2 (^ ■ Pu )(p n ■ Pl ) - m*(p, ■ Pu )} e -*&*-ft)-fa 

(18) 

It is important to notice that the kinematical regions for the momenta of 
the final state became finite for the plane waves at T — > oo due to the delta 
function which resulted from the x integration. The muon momentum of this 
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expression is that of the asymptotic state and necessarily a part of those that 
satisfy the energy and momentum conservation. Now the integrations of the 
coordinates and momenta were interchanged in the amplitude of a finite T, 
in which the states of non-conserving energy of Hq couple. In Eq. (1181) the 
coordinates are fixed and the muon momentum is not subject to the energy- 
momentum conservation. Hence the momentum is integrated over the entire 
positive-energy region as in Eq. (fTEj) . The finite-size correction is necessarily 
generated and computed from those states that do not satisfy the energy 
momentum conservation. 

Att^&e) is expressed by the four- dimensional integrals over the variable 
q = pi — p n , as the sum of the light-cone singularity, 5(Sx 2 ), and less singular 
and regular terms that are described by Bessel functions (see Appendix B). 
Thus the correlation function, A n j(5x), is expressed as 

A w ,{(fa) =2i \mlp u ■ (p w + i-^-] - 2i(p w ■ p v ) ( p n 



OCX I \ OCX 



<St) 
Air 



6(\) + I[ e9ular + I 



\ = (5x) 2 = 5t 2 - Sx\ (19) 



where j^ e9ular i s composed of Bessel functions and has a singularity of the 

form 1/A around A = and decrease as e -™v^l or oscillates as e irh ^ 1 at 
large |A| (see Appendix B), where fh 2 = m n 2 — m 2 . e(St) is a sign function 
and 5(A) is Dirac's delta function. 

Next, we substitute Eq. ffl9|) into Eq. (|T7|) and integrate over x\ and x 2 
with the center coordinate X M = (x^ + x%)/2 and the relative coordinate 
r = x\ — x*2- We have the expression from the light-cone singular term 



e(St) 



Js(x) = C 5 (A)-|^|-exp 

3 

C m = M2a \ 4> c (8t) = ujjt = ^5t. (20) 

The phase <p(5x) of the original expression became the small phase 4> c {St) 
of Eq. fTSOj) at the light cone A = 0. The next singular term is from 1/A in 



and becomes Js(x)/ \/^cr u \p u \ 2 . This is much smaller than and 
is negligible in the present parameter region. The magnitude is inversely 
proportional to \6t\ and is independent of rh 2 . This behavior is satisfied in 
general forms of the wave packets (see Appendix C). 
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The regular terms of A n j(5x) give finite contributions from the region 
f ~ 0. The first term, L\, is from j^ e9ular j n Eq. ( |19p . and behaves as 



L X = C X 



exp 



i{E v - \p u \v v )St - a u \p u \ 2 + %m\/2v v a v \p v \\St\ 



(J ( (J Tfl \ 2 3 



(21) 



in the space-like region A < 0, and L\ decreases with time as exp —mh\ \J\St\ 
in the time-like region A > 0. The second term, L 2 , is from I 2 , and is approx- 
imately the integral of exp —i — E v — a/|<? I 2 + mf^j St 

L 2 is a steeply decreasing function of \St\. 
Finally the probability is written in the form, 



m q m a range 



P = N 3 dX v dtxdt 



,i<t>c(St) 



+ 2D< h (p v )^-*(^L, 



N 3 = ig 2 ml^a u (%p n -p v /E v )V-\ 



(22) 



where Drh(Pv) and L 2 are given in Appendix B. The first term in Eq. (1221) 
oscillates extremely slowly with time St and the remaining terms oscillate 
or decrease rapidly. The integral of the first term over the finite T slowly 
approaches constant. 

1 ( d \ n 

Here we study the series } (— 2p^-p v ) n — I I L\, in Eq. (122]) . This 

/ — ' n\ \ am 2 J 

1 ( \ n 

series converges when Si = > (— 2p 7r -p u ) n — I - - - ) (m 2 )* becomes finite. 

n\ \ am 2 I 

This is satisfied in 2p n -p u < rh 2 , and the power series rapidly oscillates with 

as S2 



exp 



. Therefore the diffraction 



term exists in the region 2p n -p u < m and vanishes in outside region. 



Thus the light-cone singular term leads l-e iUJl/St and the regular term 



sv 



jreguiar ^ gj veg a g n j^ e contribution from the region r 1 ^ 0. 

The integrations over times ti and t 2 of the singular term in a finite T is 



dtidh-^-e^ 
\St\ 



T(g(u u T) - 7r) 



(23) 
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where the function g(cu u T) in the right-hand side satisfies ^(7(o; v T)|t=o = 
— u) v and g(u v T) = u u T — > oo. The last term in Eq. is can- 

celed by the integral of the short-range term i^ egular . Since g(u v T) is gen- 
erated by the superposed waves that have the form of the light-cone sin- 
gularity, and its effect remains at a macroscopic distance of the order of 
2chE u / (mlc 4 ), we call this the diffraction term. The last term I 2 in Eq. ([H?j) 
gives f a/v I dtidt 2 L 2 (St) = TG , where the constant G is computed nu- 
merically. Owing to the rapid oscillation in \8t\, the microscopic \5t\ con- 
tributes to this integral, and consequently Go is constant in T. 

Integrating over the neutrino's coordinate X y , we obtain the total volume, 
which is canceled by the factor V~ l from the normalization of the initial pion 
state. The total probability is then expressed as the sum of the normal term 
Go an d the diffraction term g(u u T): 

P = Nt J 2 ^*> BK T) + Go] . (24) 

where N4 = 8Tg 2 a u and L = cT is the length of the decay region. The 
quantity P is the neutrino flux measured by its reaction with the physical 
state of the target with finite size a v . The diffraction term corresponds to the 
energy-non-conserving term S^fT], and the normal term corresponds to the 
energy conserving term S^[T}. At T — > 00, the diffraction term vanishes 
and at a finite T, the probability contains the diffraction component, which 
is stable with respect to variation of the pion's momentum. 

3 Unusual behavior of diffraction component 

Next, we evaluate each term of Eq. (1241) . In the normal term Go, the energy 
and momentum are approximately well conserved, p n = pi +p u , and satisfies 
2p-ir ~Vv — m "i — m "i- Thus, the branching ratio of the electron mode is about 
10" 4 of that of the muon mode [221 1231 12D 123! • Integrating over the neutrino's 
angle, we find that this term is independent of a u , which is consistent with 
the condition for the stationary state [TTJ. The transition rate agrees with 
the value obtained by the ordinary method. In the diffraction component, 
g(u u T), the energy and momentum are not conserved, p n 7^ pi + p u , and the 
inner product, p-^ ■ p u , is not expressed with the masses of pion and charged 
lepton. Instead, the convergence condition requires that this term be present 
in the kinematical region, \p u \(E n — \p n \) < p^-pu < m 2 /2. Consequently the 
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Figure 2: Total detection rate of the muon neutrino at a finite distance L. 
The constant (red line) shows the normal term, and the diffraction term is 
given on top of the normal term. The horizontal axis represents the distance 
in meters, and the normal term is normalized to 1.0. The neutrino mass, 
pion energy, and neutrino energy are 1.0 [eV/c 2 ] or 0.2 [eV/c 2 ], 4 [GeV], and 
700 (blue triangles) or 800 (green circles) [MeV], respectively. The excess 
varies with the distance for m v = 1.0 [eV/c 2 ] and is almost constant for 
m v = 0.2 [eV/c 2 ]. 

electron mode is not suppressed in the diffraction component and momentum 
is integrated in this region. 

For the muon mode, / = /i, the angle of the neutrino in the diffraction 
component is slightly different from the angle of the neutrino in the normal 
mode, so it is impossible to experimentally distinguish the latter region from 
the former. Therefore, we add both terms, and the total probability thus 
obtained is presented in Fig. (F2]) for two neutrino masses m u = 1 [eV/c 2 ] and 
0.2 [eV/c 2 ], a pion energy E n — 4 [GeV], and two neutrino energies E v = 700 
and 800 [MeV]. For the wave packet size of the neutrino, we use the size of 
the nucleus having a mass number A, a v = A~3 /m, 2 . For the 16 nucleus, this 
size takes the value o v = 6.4/m 2 . From Fig. ([2]), we see that an excess of flux 
varies with the distance for distances L < 1200 [m] for m v = 1 [eV/c 2 ] and is 
almost constant for 0.2 [eV/c 2 ] and that the maximal excess is approximately 
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Figure 3: vN total cross section from MINOS and NOMAD collaborations (a) 
and total cross section of sums of normal and diffraction terms in geometries 
of MINOS and NOMAD collaborations (b). The horizontal axis gives the 
length in meters and the vertical axis gives the ratio of cross section to energy. 



The neutrino mass, 



0.2[eV/c 2 }, is used for theoretical calculation. 



20% of the normal term at L = 0. The slope at the origin L = is determined 
by u v . The diffraction term varies slowly with both distance and energy. 
For this situation, the typical length is L [m] = 2E u hc/ '(m^c 4 ) = 400 x 
E v [GeV] I ml [eV 2 /c 4 ] . The neutrino's energy is measured experimentally with 
uncertainty AE V , which is of the order 0.1 x E v . This uncertainty is 100 [MeV] 
for 1 [GeV] neutrino energy and the diffraction components of both energies 
are almost equivalent to those given in Fig. (|2]). For larger uncertainty in 
energy, the computation is easily done using Eq. (124]) . 

Fig. ([3]) compares the total cross sections of the neutrino nucleon reac- 
tions of the NOMAD [25] and MINOS J2S] collaborations. The geometry of 
each experiment is taken into account in the theoretical calculations. The 
theoretical values obtained depend slightly on energy and agree with exper- 
iments. The slight energy dependence of the total cross sections at high en- 
ergy, which is difficult to understand from the standard quark-parton model, 
becomes understandable when considering the diffraction components. 

For the electron mode, / = e, the finite-size correction becomes much 
larger than the normal term in short distance regions. This is because the 
normal term is suppressed by helicity suppression (which originates in energy- 
momentum and angular momentum conservation) [22l [23l |21~] 124"] whereas the 
diffraction term is not affected by the helicity suppression. Its magnitude 
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is about 0.2 of that of the muon mode from Fig. (J2J). Thus, the excess of 
rate becomes enormous. Fig. (Tjf]) shows the relative fraction of the electron 
mode that includes geometrical configurations of detectors. The value is 
10 -4 in the normal mode, around 0.1 at a distance of a few meters owing 
to the diffraction component, and around 0.01 at a distance of 10 to 100 
meters. The diffraction component becomes negligible in a distance above 
1000 meters. The enhancement factor is 10 3 in the former case and 10 2 in 
the latter case. For the LSND experiment [27], the distance is about 20 
to 30 meters. The theoretical values are also plotted in Fig. @ and agree 
well with the experimental values. Thus, we consider the LSND events to 
signal of neutrino diffraction. The diffraction in a configuration of the first 
high-energy neutrino experiment, TWN [SB] collaboration, is also compared 
in the same figure and an agreement is obtained. Theoretical values depend 
on the neutrino mass weakly in the range 0.2 — 1.0 [eV/c 2 ] in this figure. An 
energy spectrum at low energy region of around 100 [MeV] or less, on the 
other hand, is sensitive to the mass. Further experiments at shorter distance 
may be able to confirm the diffraction of the electron neutrino. 

Thus the probability of the detection process described by S[T] expressed 
by Eq. ( 124"]) deviates from the probability given by S[oo] and has the large 
finite-size correction. If the charged lepton is observed simultaneously, the 
detection rate of the lepton has also the large finite-size correction. This 
result is different from the probability to detect only charged lepton but 
does not contradict with these ordinary experiments because the boundary 
conditions are different in two cases. 

The finite-size correction of the probability to detect only a charged lepton 
is computed with S[T'} which satisfies the boundary condition for the charged 
lepton. Here T' = T M — T^ is the time interval for observing the charged lep- 
ton and the probability is expressed by Eq. 021]) with lo v — Y ui — mfc 4 / {2E[h). 
Since charged leptons are heavy, uiT' becomes very large and g{u{Y') becomes 

m for macroscopic T'. Thus, the probability to detect the charged lep- 
ton does not have a finite-size correction and agrees with that of the normal 
term. Although the light-cone singularity forms in both cases, the diffraction 
component becomes relevant only when the detected particle is very light. 

The probability to detect a charged lepton depends on the boundary 
condition of the neutrino. When a neutrino is detected at T„, the charged- 
lepton spectrum includes the diffraction component but, when the neutrino 
is not detected, the charged-lepton spectrum does not include the diffraction 
component. The latter result is standard, whereas the former is not, but may 
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Length between decay region and detector [m] 

Figure 4: Relative fraction of the electron mode computed theoretically for 
LSND geometry and that for TWN geometry are compared with experiments 
of LSND and TWN [28] collaborations. TWN (EXP) and LSND (EXP) give 
the experimental values and TWN (Diffraction) gives the results calculated 
with the parameters m v = 0.2 and 1.0 [eV/c 2 ], E v = 250 [MeV], and P n = 
2 [GeV/c]. LSND (Diffraction) is calculated with m v = 0.2 and 1.0 [eV/c 2 ], 
E v = 60 [MeV], and P n = 300 [MeV/c]. Flavor oscillation (T2K) shows the 
results for sin 2 6 13 = 0.11, e5m 2 3 = 2.4 x lO" 3 [eV 2 /c 4 ], and E u = 60 [MeV]. 

be verified experimentally. 

For three neutrinos of mass m Vi and a mixing matrix Ui, a , the diffraction 
term for a neutrino of flavor a is expressed as ^^(w^T)!^^! 2 , whereas the 
normal term is expressed as | ^ Ui^D(i)llj a \ 2 where i is the mass eigen- 
state, a is the flavor eigenstate, and D{i) is the free wave of m Vi . Thus, 
the diffraction term depends on the square of the average mass, m 2 , but the 
normal term depends on the mass-squared difference <5m 2 . At L — > oo, the 
diffraction term disappears and the normal term that includes the flavor os- 
cillation remains. Fig. (jl]) shows that the probability to detect the neutrino 
in the region above 1000 meters is described by the flavor oscillation and that 
below a few hundred meters is described by the neutrino diffraction. 
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We now compare neutrino diffraction with the diffraction of light passing 
through a hole. The former is that of a quantum mechanical wave and 
appears in the probability distribution. The diffraction pattern forms in a 
direction parallel to the momentum and with a phase difference u v 5t of the 
non-stationary wave. Its size is determined by u u , which is extremely small 
and stable with respect to variations in the energy and other parameters. 
Thus, without fine tuning the initial energy, the diffraction is easily observed, 
whereas the latter is that of a classical wave and appears in the intensity. 
The diffraction forms in a direction perpendicular to the momentum and 
with a phase difference cu^ B St where u^ B = cp/h for the stationary wave. Its 
shape is determined by u^ B , which is large and varies when the parameters 
are changed. Thus, the initial energy must be fine tuned to observe the 
diffraction of light. 



4 Summary and future prospects 

We found that the wave function of the pion and decay products at a finite t 
has a kinetic energy of continuous spectrum and possesses wave natures which 
are very different from that at t — oo. The rate to detect the neutrino at a 
finite-time interval T is uniquely computed with the wave packet and reveals 
a large finite-size correction. Because the wave nature causes the correction, 
that has unusual properties such as the non- conservation of the kinetic en- 
ergy, the non-suppression of the electron mode, and others. Applying S[T] 
that satisfies the boundary condition of a finite T, we obtained the results 
that can be compared with experiments. Pion and muon waves accumulat- 
ing at the speed of light form the light-cone singularity, which is real and 
extended in the large-distance region, thus strongly influence the neutrino. 
The probability to detect a neutrino produced along the light-cone singu- 
larity is subject to the finite-size correction of the diffraction phenomenon. 
The diffraction pattern is determined by the difference of angular velocities, 
oj u = uj e — ujf, B ', where u E = E u /h and uj* b = c\p v \/h. The quantity u u 
takes the extremely small value m 2 u c A / (2E u h) because of the unique neutrino 
features [SI El IB]- Consequently, the diffraction term becomes finite in the 
macroscopic spatial region of r < 27 ^"ji c and affects experiments in the mass- 
dependent manner at near-detector regions. Excess in neutrino flux in this 
region observed by K2K [29], MiniBooNE [30], and MINOS [31] may be con- 
nected with the diffraction component. In addition, the excess in electron 
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neutrinos known as the LSND anomaly may be attributed to the diffraction 
component instead of to flavor oscillation, which would resolve the contro- 
versy between LSND with others [32J. We compared all previous neutrino 
experiments and found that the new contribution in near-detector regions 
derived from the neutrino diffraction is surprisingly consistent with them. It 
would be important to confirm the neutrino diffraction with precision experi- 
ments and to find the absolute neutrino mass. If the mass satisfies ml ~ 6ml 
instead of ml 3> 5ml, the neutrino fluxes would show more complicated be- 
havior. 

We described a new quantum phenomenon for neutrinos specific to parti- 
cles of extremely small mass, and discussed the physical quantity determined 
by the absolute neutrino mass. In studying the quantum effects, we took the 
lowest order term in Gf, and ignored higher-order effects such as the pion's 
life-time, effects of electroweak gauge theory, and mean-free-path. The effects 
due to a propagator of and higher order corrections of the renormalized 
theory do not modify the amplitude in the lowest order in W's mass, Myy, 
and the light-cone singularity, hence our results of the paper are kept in- 
tact by them. They will be included along with other large-scale physical 
phenomena of low-energy neutrinos in subsequent presentations. 
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A Appendix A. The finite-size correction to 
Fermi's Golden rule 



In computing a scattering cross section and decay rate with Fermi's Golden 
rules, the following formula for a large T for a smooth function g(u), 

|^ 9M (Ef^y = T/ (i x 9 (x/T)(^) 2 = 2.T 9 (0), (25) 

is used [33, 34J. A 1/T correction to this formula is difficult to obtain with 
this expression, because an expansion, g(x/T) = g(0) + J2i fl^WAK^/T)', 
leads the integral over x for the coefficient of 1/T Z diverges for even number 
I. 

In the computation of the finite-size correction with S[T], the integral is 
regularized according to the boundary conditions of experiments expressed 
by the wave packet. Eq. (114"]) is such amplitude that satisfies the boundary 
conditions for detection of the neutrino at a finite-time interval T. Hence 
the finite-size correction is computed in a unique manner. The probability 
depends on the angular velocity oj at around a root of u — 0, where u = 
Ei + E v — E n — v v ■ (pi + p u — p w ) in Eq. f|T5|) . of the text. Due to the velocity 
dependent term in u, there exist roots of bp ^ and 5E ^ 0. 

Fig. (4) shows a solution of Sp^O. A behavior of the momentum depen- 
dent amplitude Eq. (|15p and that of the angular velocity u in the parallel 
direction to the neutrino are shown in (a) and (b), where the angle between 
the lepton I and neutrino vanishes, 8^ u = 0. The angular velocity vanishes 
and the amplitude shows the clear peak at E/H = 128 corresponding to the 
second solution, while the normal solution is not on this axis but has the 
energy E/H — 155. 

The roots for pi in the case 5p = constitute an ellipse that is almost 
symmetric of a finite size. The angular velocity u varies in the normal di- 
rection to the ellipse with a distance s linearly with a slope of the order of 
Pk or m-x- The slope is not small and the integral on the muon momentum 
along the normal direction converges rapidly owing to the coefficient T and 
the formula ^ sm(a/r/2) ^ _ 2nT5(cu) is valid. 

The roots for 5p ^ constitute also an ellipse, but the size is large and 
the shape is asymmetric. The major axis is of the order El/m 2 u and the 
minor axis is of the order E v jm v . The derivative of the angular velocity u 
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Muon Energy [MeV/ fi] 





Muon energy [MeV/fi] 



Figure 5: (a) Muon momentum dependence of the amplitude Eq. ( Tl5|) and 
(b) the angular velocity u of the text along the axis Q^ v = are given. The 
horizontal axis shows the energy of the muon and the vertical axis show the 
magnitude of the amplitude in (a) and that of the u in (b). Parameters are 
h = l i(Tv = \j m ^ C T=100 [m], E u = 30 [MeV], E w = 200 [MeV]. 

along the normal direction is computed similarly and becomes as small as 



du 
ds 



ml 
E v 



(26) 



Because this slope is so small that the integral in s does not converge fast. 
The integral at the finite T deviates from that of the infinite T and the 

2 

deviation is determined by ^T, which leads the finite size correction. In the 

2 U 

time T of ^T < 1, the deviation is sizable. 



B Appendix B. Light cone singularity. 

Innumerable states at the ultra-violet energy region in a relativistic invariant 
system lead the correlation function A n< i(5x) to be expressed by the four 
dimensional integrals of the variable q = pi — p w , which is conjugate to Sx, 
over the regions < q° and — p° < q° < 0. 
The integral over < q° is expressed as, 



d_ 

d5x 



d 
d5x 



(27) 
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where 



, 9(q^ 



4vr 4 



Im 



1 



q 2 + 2p 7T -q + m 2 — ie 



and m 



T-,2 



m 2 — mf. By expanding the integrand of 1% in p n -q, we have the 
expression with the light-cone singularity j35j|36], S(X), and less singular and 
regular terms that are described with Bessel functions, 



h = 2i 


le(6t) i 

4:71 


f short 


im 2 


8vr£ 


Dm ( -1 


a ) 

d5x 1 



-5(X) + 1\ 



regular 



2U{ 



regular 



, d 
85 x 



f, 



short? 



=eG)(-^-(^)w) , -«=*^ 



(2S 



where Ni, J±, and Ifi are Bessel functions. 

The integral over — p° < g° < 0, I2, has no singularity and is computed 
numerically. 



C Appendix C. Theorem: general wave pack- 
ets 

In general wave packets, the leading term J$/x) has the phase and the magni- 
tude of the form Eq. ( I20l) . The phase is given as a sum of 4> c ($t) and a small 
correction that becomes 0(1/ E v ) in general and 0(1/E 2 ) in time inversion 
invariant wave packets. 
(Proof.) 

For a general wave packet, the Gaussian function is replaced with the 
wave function if) (x — vt) 

tP(x-vt) = J dkdfoe'^^^^^^faik, k T ), 

dj = SyftE, (29) 

where the last term in the exponent is from an expansion of E(p + k) and 
makes the wave packet spread with time and ipk(ki, for) is a wave function 
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in the momentum. Since the coefficient CV, in the longitudinal direction is 
negligible for the neutrino, it is neglected. The probability integrated in 
center coordinate X is expressed by the square of the modules of i()(ki, kr), 
ip(ri — v u 5t, fy), in the form 

^(r, - v v 6t, r T ) = f dJfe,djfe T e ifc,(r, - ,, «' w) e < ^"^ +<( ^ /2B)ft |^ fc (A;,, k T )\ 2 . (30) 



Then the light-cone singularity leads the slowly varying probability 
J 5 (X) = J dre^W-vJtf-^-SiX) 

= vre-^ (l + 7 )^e^)^), (31) 
at a large \St\ region, where the constants are 

wo = J dk\Mh o)| 2 , 5 = § + A., 7 = *L + | (J_) 2 _ (1 _ Vv f 5t \ 

dt = wo' 1 [ dkklMkhfyl 2 , d 2 = w ~ l [ dkkf\Mk,0)\ 2 - (32) 



In high energy limit E — > oo, Eq. (1311) coincides with that of the Gaussian 
wave packet. In a wave packet of time reversal invariance, \w(ki,0)\ 2 is an 
even function of ki, hence d\ = 0. The correction terms become negligible 
for E u > 200 - 300 [MeV]. 
( Q.E.D.) 

Thus the theorem is proved. This theorem demonstrates that the prob- 
ability to detect the neutrino measured with the uncertainties of 5E ~ \5p\ 
by its collisions are modified with the diffraction component. 
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